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FOREWORD 


This report was prepared by the Polymer Branch of the Nonmetallic 
Materials Division. The work was initiated under Project 7342, ‘‘ Fundamental 
Research on Macromolecular Materials and Lubrication Phenomena,’”? Task 
No. 734203, ‘‘Fundamental Principles Determining the Behavior of Macro~ 
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Materials Laboratory, Air Force Systems Command, Wright-Patterson Air 
Force Base, Ohio, with Dr. M. T. Gehatia acting as project scientist. 
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The manuscript was released by the author in May 1969 for publication as a 
technical report. 
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ABSTRACT 


As part of a series of investigations of the fundamental properties of 
certain aromatic-heterocyclic polymers, a study has been done on the equi- 
librium sedimentation of samples of poly (2,2'-m~phenylene~5,5'~bibenzimi- 
dazole) in dimethylacetamide. Although these samples were obtained by fractional 
precipitation techniques from a whole polymer, it was found that their behavior 
was best described by assuming some samples to be composed of at least two 
major sub-fractions. In addition, this analysis required that each sub~fraction 
have a relatively strong interaction with itself and other sub~fractions, Because 
each sample consisted of a small number of major fractions a method was 
developed based upon considerations from pauci-disperse systems to charac~ 
terize molecular weights and polymer-polymer interaction terms. This report 
details the development of this method. The concentration of these samples is 
represented by: 


. 2 
c ie g,, exp [pw r-Ryo t+ (R- RO) c, |. 
where c is concentration; x = square of the distance from the center of rotation; 


h, a constant proportional to molecular weight of fraction n; R__, Ran are 


interaction coefficients and w is the angular velocity. 


The distribution of this Abstract is unlimited. 
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LIST OF SYMBOLS 
r = the distance from the center of rotation 


r me the distance from the center of rotation to the meniscus 


tr = the distance from the center of rotation to the bottom of the cell 
x =r" 

m = ra 

boos ne 

es the distance from the center of rotation to the initial boundary in 


a velocity experiment 


sedimentation constant 


D_ = diffusion constant 


angular velocity 


€ 
i] 


2 
= hw 


= time 


T = absolute temperature in °K 


=] 
tt 


universal gas constant, or interaction parameter in general 


density of solution 


D 
tt 


i = density of solvent 


Vs = partial specific volume of polymer in solution 


ec = concentration in g of polymer per g of solution 
c® = initial concentration 
H = 2RT/(1-Vp) 
M = molecular weight 
h = M/H 
hw 2. 


a constant defined for 9 ~temperature where, g=ce_ 


g = 
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LIST OF SYMBOLS (CONT) 


M. the above defined quantities corresponding to fraction n, 

n 
a 
R 


ines interaction parameter appearing in an expression for c_, and 
caused by C. : 
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SECTION I 


INTRODUCTION 


Aromatic-heterocyclic polymers are a class of interesting materials that 
are currently being developed due to their resistance to high temperature, 
Relatively little has been done toward determining many of the fundamental 
parameters that govern the physical behavior of these chains. This report 
covers one phase of an investigation of the dilute solution properties of one of the 
earlier, high molecular weight aromatic~heterocycles, poly (2,2'-m~phenylene~ 
5,5'~bibenzimidazole) (PBI), dissolved in a good solvent dimethylacetamide 
(DMAC). Specifically this report concerns equilibrium sedimentation measure~ 
ments on samples of this polymer in an effort to define molecular weight, and 
polymer-polymer interaction parameters, 


Preliminary measurements on these samples suggested that each sample 
consists of a relatively small number of rather discrete molecular weight 
ranges. In addition, it was suspected that relatively strong interactions could 
exist between various species. Therefore, early consideration was given to 
describing this situation by applying a method developed from study of pauci- 
disperse systems. 


This is discussed in the next section. 
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SECTION II 


THEORY 


A MONODISPERSE SYSTEM WITH INTERACTION 


Consider a monodisperse polymer with significant concentration dependence. 
If higher than first order terms can be neglected, the concentration achieved in 
equilibrium-sedimentation can be expressed by the following equation: 


hw*c- Re 


c = ge (1) 


The relationship expressed by Equation 1 and especially the interaction param~ 
eter R were discussed by Fujita (Reference 1), Casassa (Reference 2), Gehatia 
and Wiff (Reference 3), and others. 

By differentiating Equation 1 with respect to x one obtains: 


Se = (nw?-r $2) c (2) 


which leads to the following working formulas: 


-lidc . 2_ _dc 

c aK hw aa (3) 
| es 2 dc -t = 
c hw (<= ) R (4) 


A plot of c = a vs. se , according to Equation 3, should give a ee line with 


a slope equal to (-R). aaa ei according to Equation 4, a plot of c_ * vs. ( ey ag 
gives a slope equal to hig” » Which is a quantity proportional to the molecular 


weight (References 3 and 4). 


By knowing h Pe and R one can evaluate the constant g from Equation 1 and 
thereby fully characterize the system under consideration. 


SS Sa saan 
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SYSTEM OF TWO DISTINCT INTERACTING FRACTIONS 


Consider a solution of polymer comprised of two interacting fractions (1 and 


2). The concentration can be expressed as: 


c = Cc, + Co (5) 
where: 
2 
c = ge" x — Ry 46, - Ryace (6) 
{ 
and: 2 


It has been assumed that in acertainregion (m £ x S x*) the fraction 1 prevails, 


e.g., that c, ~ candc, << c,. Since c, =c-c, the exponent in Equation 6 can 
1 2 1 1 2 


be modified: 


“Re. >R c, = -R pe ~ CR, 


22 , a R, de (8) 


F 2 


and Equation 5 can be approximated by the formula: 


2 "i 
~h,wx-R 
c® ge ! moe [ 1-(R, 2- Ry) ea] + Cy (9) 
Denote: 
2 
h,wx - 
ge OR 2g (10) 
and: 
Rie RF Kn . (11) 


Equation 9 will now lead to the following expression for Co} 


c- G, 


I-KiG, 2 


The approximation is justified if the quantity Knte << 1, e.g., if, for a 
certain x> x*, Cy is not negligible in comparison tobe c, but is still small 


enough to make the following transformation valid: 
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For a known value of Kn one can evaluate the c 
Kine << 1. According to Fujita (Reference 1), 


curve in a region in which 


2 


R = Mf (14) 
and: 
R =M, f (15) 


where fie and toy are the cross coefficients of interaction and: 


fi,eF fo, (16) 
Using this assumption of a symmetric interaction one obtains: 


Mo h w* 
R, = R +) =R 2 (17) 
2,! t,2 (a2) ek no?) 
By inserting this expression for Ry 1 into Equation 7 one obtains the following 
formula: 


a ; how how 
Inc, = Ing, - 4, R c +[ 2°) R —R Je (18) 
h, a t,2 h,w 1,2 2,25 2 


Equation 18 is an expression for Cy with three unknown parameters; Ing,» 
2 2 _ 2 2 = 
Oyo *Pagu®) Rg =a, and [ gw™/oyu"), 9-Ry p] =? 


For a large number of c values one can determine the parameters in 
Equation 18 as well as the total error of the system. In this calculation a 
measure of the error was taken as: 


2. 2 
4, * » oni (19) 
where: 
8n,i° In (Calm im In(galyt a_(c)it bolce Mn,i (20) 


One must remember that such a treatment was originally suggested for 
the case when Kn is a known quantity. In a real case such a value is not known, 
a priori, However, the following procedure can be applied. A set of values can 
be assumed for Kn and accordingly A an canbe evaluated, The ‘‘best fit,’’ e.g., 
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the minimum of a plot A a vs. K_» gives an acceptable value for K_,° In such 
a case the values of Kn as well as other parameters, are determined. These 
parameters finally lead to evaluation of the following constants: M,> M,,; gy 
Eos Ry iy Ry 2? Roy and Ro 2 which fully describe the system comprised of 
two interacting fractions (Figures 16 and 17), 


The search for appropriate Kn values can be accomplished by use of a 
high speed computer. In case of a few fractions one can find more than one 
distinct Kn value, i.e., {k,, ; with the corresponding set M,» g. Ri R a’ 
and Ran If there are too many fractions in the system Kcn values cannot be 

eet ? 
distinguished by searching, and ‘‘noise’’ is created. 


In case of two fractions, and if allparameters of the transcendental equation, 
Equation 7, are known, ¢, can be precisely evaluated for all values of x without 


need of an approximation. 
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SECTION III 


EXPERIMENTAL 


Bulk PBI was dissolved in DMAC and fractionated with hexane. The fractions 
were purified in the following manner, Each sample was redissolved in DMAC, 
precipitated and washed with methyl ethyl ketone (MEK), filtered and washed 
again with MEK. The process of washing was continued; however, MEK was 
replaced by mixtures of MEK-MEOH with successively decreasing amounts of 
MEK. The polymer was further washed with pure methanol and replaced by 
mixtures of MEOH~ether with successively decreasing amounts of MEOH. 
Finally, the remaining polymer was washed with ether and dried under vaccum 
(this method of purification and fractionation was suggested by T. E. Helminiak), 


A set of sedimentation equilibrium experiments was carried out with four 
samples of PBI in DMAC at 40°C. An aluminum, 12mm, 4°, single sector cell 
containing the solution and another similar cell with DMAC were inserted 
into a J-rotor. The resulting Schlieren curves appeared, therefore, with a 
base line in addition to the sedimentation curve (Figures 2-8). 


Each sample was measured at different rotor speeds. It required about 
14-17 days to achieve the first equilibrium. If the speed was decreased the next 
equilibrium could be achieved within 7-10 days. However, it took only a few 
days to achieve equilibrium when the speed was increased rather than decreased. 
This surprising behavior, which contradicts theoretical considerations and ex~ 
pectations, has not yet been explained. In several cases no equilibrium was 
achieved as the speed was lowered and the Schlieren curve dissipated, 


The quantity cs has been determined from auxiliary velocity runs carried 
out with a synthetic boundary cell. 


Finally the plates were enlarged and the coordinates a and Ar were 


measured, 
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SECTION IV 


COMPUTATION 


The experimental data was analyzed by applying Equations 3 and 4, 
(Figures 2-16 and Tables I-X]). As one canreadily see no linear plots have been 
obtained, Therefore, the samples under consideration are not homogenous. On 
the other hand, one can also observe that these plots clearly show a straight 
line in a zone close to the meniscus (m< x < x*, where x* is a special value 
different for each sample and speed). Only sample three led to a straight line 


over all values of x, 


Such a peculiarity of the plot can be explained by assuming the existence 
of a distinct low molecular weight fraction. This, as well as other considerations, 
suggests that there exists a very pauci-disperse system, i.e., that each sample 
is comprised of a few fractions (as a matter of fact, 2, 3, or 4 fractions). In 
addition, the low molecular weight fraction is very distinct in a certain zone 
near the meniscus, Crotal © °1 and Cc, #1 << Cy. 


The results of applying a method based upon these considerations to the 
four samples is summarized in Table XIL An equilibrium sedimentation ex- 
periment was made at each speed indicated, The first fraction parameters 
(g)> M,, and Ry) were evaluated by the procedure outlined above. The other 
parameters Gn My Rim’ and R, re were determined by finding minima as 
indicated in Figure 22, All minima corresponding to fractions within a sample 
are indicated in Table XIL 
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SECTION V 


DISCUSSION OF RESULTS 


The computation just described cannot be considered as completed. In the 
case of four fractions: 
= 2 
Cy = Gn, exp (how -Ric-R,oc,-R 


n3°s Rng, ) (21) 


each curve can be evaluated within a zone, where Kincm <<1, and the corre- 
sponding parameters can be determined within such a zone. 


However, an assumption has been made, according to Fujita (Reference 1), 
that: 


R . R 
ok = ko (22) 

Mp Mk 
which is not readily apparent (Reference 2). Without using the relationship de- 
scribed in Equation 22, the cross-~-coefficients Run (k #1) cannot be easily 
evaluated. 


The parameters hw 2 and Ri were determined from the linear portion 
of original plots. However, such an evaluation may introduce a numerical error 
which exceeds the tolerances required by the Kan’ A mon analysis, Therefore, 
a variance in the values of h,w and Ria has to be taken into account to 


1 
better fit the system. 


Finally, the analysis has been performed by comparing the lowest fraction 
of a given sample with the other fractions in the sample. The results of applying 
this analysis to four samples is given in Table XIL Samples 1, 2, and 4 show 
that they are composed of about four sub-fractions. These sub-fractions have 
molecular weights of approximately 1,500; 30,000; 60,000; and 120,000. One 
sample, namely number 3, appears to have a very narrow distribution. Its 
weight average molecular weight is about 4,500. 


ae 
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It is also significant that the diagonal elements im) of the ‘‘interaction’’ 
matrix were all found to be negative. Only values for Rim are given in Table XII, 
since no proper method was applied to determine the off~diagonal elements R nk 
for n #1, 


It has been suggested that some of the molecular weights indicated 
(Table XII) by this experimental analysis are higher than would be expected 
from condensation polymerization, Therefore, further investigations must be 
initiated for better examination of the system. Also, additional work will be 
required to more fully evaluate the validity of this method and to make neces~ 
sary modifications to include in the calculation of all R. nk parameters. 
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SECTION VI 


COMPUTER PROGRAM 


COMPUTER SYMBOL IDENTIFICATION 


Following is a list of symbols used in the computer program, Where pos- 
sible the symbols are identified with the previously derived theory. One must 
remember that the computer program was written only for the case of two 
major distinct molecular weight components. 


RECQ@RD - Used as “flag’’ for subroutine calling order. 
PNAME1 and PNAME2 ~ MWentification of experiment. 
NMAX - Number of data points read into machine. 

CO ~ Initial concentration of experimental solution, 


R1 - First radial value, The distance from the center of rotation in the 
meniscus. 


DLR - Incremented value for radial distances, Radial distance to the 
bottom of cell is given by Rl + (NMAX~1)*DLR. 


R(D - Array of radial values 
X(D) - Array R(I)**2 


CR()) - Array containing experimentally measured ordinates from 
Schlieren Curves. Once these are read they are multiplied by 


SCALE to obtain the true de, dr values. 
CX() - Array (/2r) (°/ar) 
x 
C(I) - Array of concentrations obtained from a dx=c- Cue Then 
b m 
use is made of the equation f cdx = c°(b-m) to calculate Cw? the 
m 
concentration at the meniscus, 
ZLC(D - Array 1n (C(D) 
CX@VX(I) - Array (dc/dx)/C 


XCX(D - Array (dc/dx)} 


AFML~TR-69-235 
Part I 


XC(D - Array (cv) 
NMAX1 - First point at which Cy is small but must be considered. 


NMAX2 - Highest point at which c, is small but not zero. 


DLK ~ Ri 1 + R, 2 
DKMAX - Maximum value of DLK, 


2 


XDLK - Increment value for DLK. 
2y ~ 

G() ~ Array {ge og Bl \ 

HIGMSQ ~ hw 

Rll - Ri 

Gl- 8 


Y(L,J) ~ Array containing three consecutive sets of values (three 
successive DLK values) for YN(]. 


YN - (CQ) - G@)/(-DLK*G@) 
U@ ~ Array (K - (Ryp/h,w *)*c) 


Iy(1,J) ~ Array containing the values of NMAX1 and NMAX2 corre~ 
sponding to the array Y(I,J). 


Z() - Array In Co) 


R22 - Roo 
R12 - R, . = DLK- Ry 
G2 = go 


H20MSQ - how? 
XLNG2 ~ In 5) 
XN - Number of data points sampled (NMAX2-NMAX1 + 1) 


| NMAX2 R 
DELTA = aN > { in (x) -[Ing + hw (x, - He c(x.)) 
i=NMAX1 on oe ‘ hw! 
Ri,2 


The flow of computer data is shown in Figure 1. 
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COMPUTER PRINTOUT 


S$IBFTC MAINF DECK 


fy 


3000 


7000 


COMMON X9CXsCsNMAX 

DIMENSIONS FOR COMMON 

DIMENSION X(100@) »CX(100) 2C(100) 
DATA PlsP29P3/6HPART 1s6HENDRUN 96H 
DATA P4&/6HPART 2/ 
READ(59101) RECORD» PNAMEI sPNAME2 
FORMAT (A624X 92A6) 
IF(RECORD -— P1) 2591092 

IF(RECORD -— P2) 3398993 

IF(RECORD - P4) 491194 

IF(RECORD - P3) 9491999 | 

WRITE (631000) 

FORMAT (1H1/1HA) 

WRITE (692000) PNAME1 »PNAME2 
FORMAT(1H 95CXs25HIDENTIFICATION NUMBER IS »92A6) 


‘CALL FRCINT 


WRITE(691001) 

FORMAT (1HA/1HA) 

GO TO 1 

WRITE (691000) 

WRITE (632000) PNAME1 sPNAME2 

CALL CTWO 

WRITE (691001) 

GO TO 1 

WRITE (633000) RECORD 

FORMAT(1H »36HWHAT DO WE DO WITH THE CARD LABELED »A6) 
GQ TO 1 

WRITE (631000) 

WRITE(6s7000) 

FORMAT(8(11H END OF RUNs3X)/1H1) 
STOP 


END 
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COMPUTER PRINTOUT (CONT) 


$IBFTC FRCIN DECK 


191 
109 


192 


1 


2 


2909 


3 


2001 


4 


2002 


SUBROUTINE FRCINT 

COMMON X9CX9C ¥NMAX 

DIMENSIONS FOR COMMON 

DIMENSIUN X(100) 9CX(100) 9C (100) 

DIMENSION CR(100) sR (100) sZLE (100) sCXOVL( 100) 9XLX( 100) 9X01 100) 
READ(59101) NMAXsCOUsSCALEsR1IsOLR 

FORMAT (1391P2E116420PF1005s0PE 1006) 

FORMAT(AS) 

READ(S3192) (CR(I)soI = 1LseNMAX) 

FORMAT(16F5.0) 

NO 1 1 = l»yNMAX 

R(1) Rl + FLOAT(I-1L)*DLR 

X(1) R(1)#*z 

CR(I) = SCALEXCR(1) 

CXC1) = CRED) (20%*R(1)) 

CONTINUE 
A = OO. 
COEF = Oe 
NA 2 JT = LoNMAX 

TF (Te¢E&Oe1) GO TO 2 

DLX = X(T) — X¢I-1) 

AVGCX = (CX(1)+CX(1-1))/26 

COEF = COEF + DLX*AVGCX 

AC = (COEF + C(1-1))/26 

A = A + DLX¥AC 

C(I) = COEF 

DIFsmM = X(NMAX) - X21) 

CM = (DIFBM*¥CO — A)/DIFBM 

WEITE (692000) COsDIFKMsCM 

FORMAT(1H »933HTHE CONCENTRATION FOR THIS RUN = 96€11¢4/43H DIFFEREN 
1CE BTWN SQS OF BTM AND MENISCUS = 9£15¢8/21H CONC. AT MENISCUS = 94 
2F11.4////) 

DO 3 T = LeNMAX 

C(1) = CCl) + CM 

ZLC(I) = Os 

IF( COL) oLEeO6.) GU TO 3 

ZLC(I) = ALOG(C(1)) 

CXOVC(T) = CXTI)Y/C(T) 

XCX(1) = Lle/CX(1) 

XC(I) = 1e/CC(1) 

CONTINUF 

WRITE (692001) 

FORMAT(LH s1XslHIsSXs1HRellXslLHXsl1lXs2HCRelCXsZHCX slOXs LHC s9XsSHEN 
1(C) 9 7X 94HOX/C 9 BX 9 4ri lL /CX 9 BX 9 BHI/SC//) 

DO 4 I = l»sNMAX 

wRITE (692062) IsR(1) 9X1) oCRUT) sCX(1) oC (1) sZL CCL) sCKXUVE(1) »sXCXU 1) 5 
UXCCT) 

FORMAT(LH 9 1295X9S(1PEl1e¢491X) ) 

RETURN 

END 
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COMPUTER PRINTOUT (CONT) 


“E€LTBFTC CTAU DECK 
SUBROUTINE CTWO 
COMMON X»CXsC»sNMAX 
C DIMENSIONS FOR COMMON 
DIMENSION X(100) eCX(100) sC(100) 
NIMENSTON G(109) sYN(1G0) 9li(100)sZ(100) 
READ (53100) NMAX1sNMAX2 s9KMAXsH1OMSWsR1L9G1»DKMAX sXDLK 
INO FORMAT (312393F10.39F150¢8sF1005) 
~ DO 1 T = NMAX] sNMAX 
AL = HIOMSQ*X(1T) - R11*¥C(1) 
A2 = EXP(AL) 
1 G(1) = G1*A2 

DO 2 K = 1»sKMAX 
K1 = NMAX1 
K2 = NMAX2 
NLK = NKMAX + XNLK#FLOAT(K=1) 
WRITE(6s1901) OLK 

1201 FORMAT (IH »89OXsSHKN = 9F1Q65) 
P12 = R11 - ALK 

98 CONTINUE 

WRITE (622000) NMAX1 »NMAX2 

2009 FORMAT(1H s8HNMAX1 = s12s9H NMAX2 = 612) 


DO 2 I = NMAX1l»NMAX2 
Al = CCl) = G1) 
A2 = le + DLK*G(1) 


YN(T) = AL/SA2 
TFCITSEGeNMAX1) GO TO 20 
EN TO 2] 
2h AS = X(1) - R12#C(1)/H1OMSO 
XN] = AR = JT, 
21 CONTINUE 
MCT) = X(T) = CRI2*CCI)/HIOMSQ) - XNI 
XXN = le 
LFCYN(C i) eLT*eOeLE-05) GU TU Gy 
TF CT eEQeMNMAX1) GO TO 22 
GO TO 23 
22 A2 = ALOG(YN(1)) 
XN2 = ABS(A3) ~ 1. 
22 CONTINUE 
Z(1) = ALOG(YN(T)) + XN2 


2 CONTINUE 


fi = Ne 

A2 = Oe 

AR = O, 

A4 = Qe 

A5 = Ce 

A6& = O« 

Ay = Ne 

AR = Of 

XN = FLOAT(NMAX2 -— NMAX1 + 1) 
DO 4&4 T = NMAXI se NMAX2? 

Al = Al + Co1¥U(T) 

AP = AP + (O41*80U01) )*¥¥2 
AZ = AP + LOOW¥YNCT) 
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COMPUTER PRINTOUT (CONT) 


AG = AG + (1006#YN(1j))*¥2 
AS = AS + Z(T) 
M6 = £6 + 10e*U(T)*YN(T) 
AT = AT + O61#U(1)¥*Z(1) 

4 AB = AB + 100e*YN(1)*Z(1) 
R] = A2¥AA — AB¥AG 
R2 = Al¥A4 - ABA 
(AD = AQHAD = ALHAG 


DENOM = XN*B1 -— Al*B2 — A3*83 
ANUM = A5*Bl -— A7¥B2 - AB*S3 
= AT¥AG — ABXAG 
R2 = AS*AG — AZ¥NB 
Q2 = ARHAT — ARXAF 
BNUM = XN¥®B1l = AL¥B2 — A3*53 
2] = AT¥AG — A2#AS 
R2 = AS¥AG — A1#AB 
R2 = ALHAT — AS#AD 
CNUM = XN¥B1 = A1¥*B2 — A3*63 
H2OMS@ = O+61*(3NUM/DENOM) 
B = 100¢#(CNUM/DENOM) 
R22 = 3 + R12*HZOMSO/H10MSQ 
XK = CNUM/DENOM 
XLNG2 = XK — XN1*¥HZOMSQ - XN2 
WRITE(693001) DENOMsANUM sBNUM» CNUMsK22 sHZ20MSu 9 XLiNGe 
2001 FORMAT(IH 97(1XsE156891X)) 
TF (A2S(XLNG2).GTe886) GO TO 40 
G2 = FXP(XLNG2) 


49 G2 = (Ce 
41 CONTINUE 
COEF = Ce 
DO 5 I = NMAX1 »NMAXzZ 
Pl = XK 
R2 = H20MSG#U(T) 
O34 = RXYN(T) 
5 COEF = COEF + (2(1)-81-32+B3)**2 
DELTA = COEF/XN 
TFC HZOMSQeLTAN0e) GO TO 32 
WRITE (691900) KsR129G2sH20MSo9R22sDELTA 
1000 FORMAT(1H 95XsI297H R12 = slPE10e396H GZ = s1PE106¢3:10H HeUnsa = 5 
LIPEiNe397H R22 = »91PE100399H DELTA = »o1PE1063/) 
32 CONTINUE 


NfAXL = KI 
NMAX2 = K2 
GC TO 2 


99 CONTINUE 
NMAX] = MYAX] + 1 
TEC CNMAX2-NMAX1) 6GTe2) GO TO 97 
N“AX? = NMAX2 + 1 
TF (NMAX2¢EGeNMAX) GO TO 88 
97 GO TO 98 
2 CONTINUE 
88 RETURN 
END 
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TABLE XII 
FINAL RESULTS OF THE ELEVEN EXPERIMENTS REPORTED 
ROTOR SPEED . 

SAMPLE (RPM) Mx107> gy “Raz Mylo g, RL, 
1 10, 589 1.4 3.393x1072 6.16 = 28.9 -2.135x107” 134 
50.2 715x10719 165 
97.7. 9.757x10727 89 
1 8,766 1.4 4.015x1072 6.24 30.1 1.022«107> 386 
38.6 8.520x10~6 325 
67.4 1427x1078 208 
82.6 4,.047x107 166 
105.5 5.114x1072 110 
112.9 2.545x107 95 
1 7,477 1.7 4.430x10"2 5.90 29.6 1.413x1074 241 
59.3 1 . 301x107? 217 
131.3 6.394x107 127 
2 17,250 1.6 3.520x1072 20.0 5.2 8.288x107> 75 
43.8 1.187x10729 — 186 
116.3. 6.127x107%8 ~=—_255 
2 10,589 1.6 1,050x1072 20.0 44.5 1.381x1077 122 
58.6  1.222x10™ 45 
62.1 6.474x1079 27 
66.6 2.874x10~9 6 
2 7,447 3.3 1.503x1072 13.2 32.7 5.616x107> 215 
114.3 2.545x10-8 86 

3 10,589 44 1.125x10°2, 5.3 

3 8,766" 4.7 1.966x10-2 | 4.7 No other significant components 

found for this fraction. 

3 7,447 4.2 2.740x10-2, 5.5 
4 13,410 1.6 1.901x10-2. 7.0 6.2 1.729x1074 164 
10.7 2.003x107> 160 
139.3 1. 580x1073? 117 
4 10, 589 1.6:  2.422x1072 8.3 3.2 3.352x1073 136 
9.8 2.337x10™ 38 
44.6  1.525x10-6 373 
123.1 1.190x107!9 1,558 


The subscript m(m = 2,3,...) indicates other fractions found. 
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Figure 2, Schlieren Curves for Sample 1 at Rotor Speeds: a) 10,589 RPM; 
b) 8,766 RPM; and c) 7,447 RPM 
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Figure 3, Schlieren 
b) 7,447 
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Figure 8. Schlieren Curve for Sample 4 at Rotor Speed 10,589 RPM 
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Figure 17, Error hoor Vs, (-K,) for Sample 2 at 17,250 RPM, Indicating a 
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Figure 18. Example of How More Than One Fraction is Predicted (Each 
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